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I. INTRODUCTION 

Experimental realization of Bosc-Einstcin condensates 
(BEC) in ultracold trapped gases opened up a rapidly ex- 
panding field of studies of quantum-degenerate systems 
HIj 0, U- Among the others, statistical properties of the 
condensate, such as its fluctuations and correlations, have 
attracted a wide attention. While the theory on this sub- 
ject is well developed (see e.g. Q and references therein), 
there are only few experiments that address this issue. To 
date, only the fluctuations of the total number of atoms 
in a condensed gas have been measured Q , and the sub- 
poissonian scaling has been observed. The second or- 
der correlation functions, that are directly connected to 
the condensate atom number fluctuations, have been in- 
vestigated experimentally in the collisions of metastable 
helium condensates Q and in the expanding rubidium 
condensate 0. 

One of the potential tools to measure the statistics of 
quantum-degenerate gases is based on atom-light inter- 
actions. This possibility has been noticed already some 
time ago, and has been proposed for a detection of the 
Bose condensed phase [1, |^ O, [HI , superfluidity in Fermi 
gases [H, [13, [3 [l^ and, quite recently, for a detection 
of quan tuin phases in ultracold gases in optical lattices 
[l^llTl [Tsl . [l9j . The optical imaging techniques have al- 
ready been used to measure coherence properties of BEC 
in the Raman superradiant scattering [20| and the second 
order correlation functions . 

The light scattered from a quantum gas carries infor- 
mation on atoms statistics, and thus can be used to mea- 
sure the condensate fluctuations [2l| . In the case of BEC 
in a trap the profile of the scattered light is dominated 
by a component, which depends on the mean occupation 
number of the condensate. In order to detect a much 
weaker component resulting from fluctuations, one has 
to resort to the variance of the number of scattered pho- 
tons, which can be difficult to measure. 

The situation changes, however, in the presence of a 
periodic potential. In this case the dominating classi- 
cal component exhibits interference pattern characteris- 
tic for the Bragg scattering, and the quantum component 



can be measured at the angles corresponding to destruc- 
tive interference, where the large classical component 
vanishes. This property has been first noticed by Mekhov 
et al. [Ill , and these authors have proposed a method of 
probing the statistics of an ultracold gas in a lattice , 
based on the relatively strong coupling between atoms 
and light modes of a cavity. In this case one should be 
able to perform non-demolition measurement allowing to 
distinguish between superfiuid and Mott-insulator (MI) 
quantum phases at temperature T — 0. 

In this paper we study a less complex situation of 
measuring a quantum gas statistics based on the far-off- 
resonance light scattering from a Bose gas in a lattice, 
focusing on the effects of statistics at finite temperatures. 
In order to avoid atom losses and suppress a possibility 
of perturbing the quantum state by the probing laser, 
we assume that the probing light is sufficiently weak and 
far-detuned. We show that the mean number of photons 
detected at some special angles, carries enough informa- 
tion not only to distinguish between different thermo- 
dynamic phases of the gas but also to directly measure 
the effects of the on-site atom statistics driven by quan- 
tum and thermal fluctuations. Hence, it allows to verify 
the validity of some well-grounded literature approaches, 
such as the Bogoliubov method, to describe higher order 
correlation functions in an interacting Bose gas. 

Our paper is organized as follows. In section |TT] we 
develop a model for scattering of light from ultracold 
atoms, showing that the number of scattered photons 
is directly related to the second-order correlation func- 
tion. In section IIIII we tailor our model to the external 
potential created by an optical lattice. The scattering 
from atoms in one dimensional (ID) lattice is considered 
in section IIVI where for simplicity we focus only on the 
zero-temperature statistics discussing the effects of differ- 
ent approximations. The finite temperature statistics of 
a Bose gas in a lattice is analyzed in section|Vl Sectioir lVll 
investigates scattering from atoms in three-dimensional 
(3D) optical lattice at finite temperatures. We conclude 
in section IVIIl and finally two appendixes present tech- 
nical details related to the influence of non-local Franck- 
Condon coefficients (Appendix A) and optimal configu- 
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ration of a probing laser and a photon detector in the 3D 
lattice case (Appendix B). 



II. INTERACTION OF LIGHT WITH MANY 
ATOMS 

In this section we consider a general problem of light 
scattering from a gas of bosons in an arbitrary external 
potential. We assume that the trapped atoms are illu- 
minated with a weak, and far-detuned laser light. The 
angularly resolved scattered light is measured by detec- 
tors in the far-field region. The full Hamiltonian of the 
system consists of the following parts: 



Ti ~Ha+ Ti-f + Hal + Haf 



(1) 



where Ha is the atomic Hamiltonian, Hf represents vac- 
uum modes of the electromagnetic field (EM), Hai de- 
scribes interaction of atoms with the laser light and Haf 
interaction of atoms with vacuum modes. 

The atomic Hamiltonian can be split into two parts 



Ha = Ho + Hi 



(2) 



where Ho describes the system of two-level atoms in the 
second-quantization formalism (23| 
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and the part including the atom-atom interactions reads 



Hi 



nmpq 5n5rn.9pffq 



(4) 
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Here, (gj^) is the annihilation (creation) operator of 
an atom in the ground electronic state and state n of the 
center-of-mass (COM) motion, and em (efn) is the anni- 
hilation (creation) operator of an atom in an electronic 
excited state and state m of COM motion. The oper- 
ators obey the standard bosonic commutation relations: 
bn,3m] = and [en,ej:„] = J„.m- The corresponding 
eigencnergies of atom COM motion are denoted by ?iaj^ 
and tvjj'^ for atoms in the ground and excited electronic 
states, respectively. The matrix elements ?7nmpq of the 
interaction Hamiltonian read 



dM*r.{rW^{r)M^)M^) (5) 



where we model short-range interactions through a con- 
tact potential with s-wavc scattering length as and a 
mass of the atom m. We neglect ground-excited and 
excited-excited atom interactions, assuming that for a 
weak and far-detuned probing light excited atoms con- 
stitute only a small fraction of the whole sample. 

The Hamiltonian of the EM field takes a standard 
form: 



Hf = ^ j (Pk fiwkalA^kA 



(6) 



with flkA (flkA) being an annihilation (creation) operator 
of a photon with a wave vector k and a polarization A. 

The interaction of atoms with a laser beam is described 
as follows: 



Hal = ■^^(n,5|ukjr)|m,e)e'"^*fftem + /i.c. (7) 



where we treat the macroscopically occupied laser mode 
classically. Here, Wk^ (r) characterizes a laser mode with a 
wave vector k^, , uj^ is the laser frequency, f2 is a Rabi fre- 
quency of the atomic transition, and the Franck-Condon 
coefficients (n, (^juk^ (r)|m, e) describe a transition am- 
plitude between COM motion states n and m of the 
atoms in the ground and excited electronic states, re- 
spectively. Typically, for a single probing laser, we have 
WkL (i") = e^^^^ (running wave), and for the two counter- 
propagating probing beams Uk^ (r) = cos(kLr) (standing 
wave) . In general uv, (r) can also represent the modes of 
an optical cavity [l6j . 

The part of the Hamiltonian that describes coupling of 
atoms with quantized EM field is given by: 



Haf 
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(8) 



in which CkA = i/wk/ (167r'^eo?i) (d • CkA), d is a dipole 
moment of the atomic transition, Mk(r) is a mode func- 
tion of the EM field with a wavevcctor k and frequency 
(jJk, and ekA is a unit vector perpendicular to k describing 
the mode of light with a polarization A. 

We solve the quantum equations of motion in the 
Heisenberg picture under the following approximation: 
i) we assume that the atomic operators arc driven only 
by the dominating laser mode of the EM field, neglecting 
the back action of atoms on the laser mode, ii) the quan- 
tum dynamics of the vacuum modes is determined by the 
evolution of atomic operators, ignoring the back-action of 
the vacuum modes, which is equivalent to neglecting the 
process of spontaneous emission, iii) for the weak and 
far-detuned laser field we perform adiabatic elimination 
of the weakly populated excited state. Our approxima- 



tions are analogous to those used in [21|, with the only 
difference that here we perform the adiabatic elimina- 
tion of the excited state, instead of assuming short prob- 
ing pulses. We carry out our derivation neglecting the 
interactions between atoms and we comment on the gen- 
eralization to the interacting gas case at the end of this 
section. 

The equations of motion for the atomic operators in 
the interaction picture with respect to Ho'- 5m (i) = 
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5m(i)e" 
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and en{t) = en(i)e-*('^" read 
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.9m (t), 
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where ??nm(k) = (n, e|uk(r)|ni, (7), A = — and we 
have introduced rcscalcd time variable r = At. We solve 
Eqs. ([U and pO|) by applying the Laplace transforma- 
tion F'^(s)= dr e^'''^ F{t). The Laplace transformed 
equations take the form 



■S5m(s) - 5m(0) = 

-«^5I^nm(kL)Sn(s) 



(11) 



^(s)~e„(0) = 
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For a far-detuned light the prefactor on the right-hand- 
side of the equations is small: r2/A <C 1, and, in princi- 
ple, the equations can be solved by iterations in a pertur- 
bative manner. Here, however, we proceed with solving 
Eq. (fT2|) for e^(s) and then substituting the result into 
Eq. (fTT|) . For a far-detuned light we apply A ^ uj^,uj^ 
and we use the identity ^ri^^(kL)rinm'O^L) = ^mm', 

n 

which results in 



We substitute Eqs. (|14p and ([15)) into equation of mo- 
tion of the E-M field operators akA(i) = akA(06~''^''* in 
the interaction picture. In the lowest order in 51/ A this 
yields 

di^xit) - akA(O) = 

nn'm 

X J sm 1 (16) 

where oj^n' = + w^, — (jJ^- At t ^ 00 the sine term 
will produce a term proportional to the delta function, 
describing the energy conservation in the process of a 
single photon scattering: Wk = Wk^ + uj^, — w^. How- 
ever, in our case we are interested in the total number 
of photons scattered into a given solid angle, and not in 
the spectrum of the scattered light. Hence, we use the 
approximation oj^n' ~ ^l- This condition is also appli- 
cable in the physical systems where the natural linewidth 
r associated with the atomic transition is broader than 
frequencies of atom COM motion: F ^ u;^. 



Now, by using Eq. (|16p and approximation Lo^n' ~ 
we calculate the mean number of photons with a wavevec- 
tor k and a polarization A 



ftkA (0 OkA [t) 



(17) 



where the function F(k, k^) is defined as follows: 



F{\.M)^Y. (n|u^(r)^kjr)|n')(m|^k(r)u^jr)|m') 



s -\- 1 



-l=p ( 5m(0) - Z— 2^77nm(kL)— - 



(13) 

By substituting back this result into Eq. ([T^ . and per- 
forming the inverse Laplace transformation we obtain the 
following time-dependence of the atomic operators 



5m (i) = 5m(0)e 



(14) 



e„(t) = e~„(0) 
+ 2A 



^'7nm(kL)gm(0) 
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Here, ljac = denotes AC Stark shift of atomic levels 
in the field of the probing laser. In Eq. we have not 
included terms of the order of fl/A , which are propor- 
tional to en(0), since they do not give any contribution 
to the mean number of photons, assuming that there are 
no excited atoms at the beginning. 



X (gU0)5n'(0).9L(0)5m'(0)) 



(18) 



Notice that, in Eq. pS]) all the matrix elements are calcu- 
lated between COM states of ground-state atoms |n, g) 
and to shorten the notation |n) = |n, g) . In the particular 
case when the mode functions Mk(r) and Uki(r) are the 
plane waves, i^(k,kL) reduces to the Fourier transform 
of the second-order correlation function in atomic field 
operators *I'g(x) of the atoms in the electronic groimd 
state 



F(q) = Id^x 



(19) 



where q = k — k^ is the wave vector of the momentum 
transfer. In the rest of the paper we will use the F{q) 
function only. 

An analogous result is obtained when considering the 
scattering of neutrons from liquid helium (23 |. In that 
case the number of scattered particles associated with 
the momentum transfer q and the energy transfer to the 
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system hu; is described by the dynamic structure factor 

X (*B|/5(x),5(iJ -E- M/5(x)|^'b) (20) 

where p(x) = *I'^^(x)^(x), H is the Hamiltonian of the 
system, and I^'b) is an eigenstate with energy E. By 
integrating over energies of the scattered particles one 
obtains the static structure factor 



S{q) 



duj S'(q, Lo) 



(21) 



which is equivalent to our function -F(q) describing an 
amplitude of scattered photons integrated over photon 
frequencies [l^. We will refer to -F'(q) as the structure 
function in the rest of the paper. 

For evolution time t much longer than the time scale 
determined by the optical frequencies wl, we can apply 
the following identity 



lim 

t — >oo 



sin^ ((tJk — wl) t/2) TT 



(22) 



to show that for the weak and far-detuned laser the scat- 
tered light described by Eq. ([T7]) has spectrum centered 
around elastic component. 

The total number of photons scattered into a solid an- 
gle dil is equal to 



dN, 



ph 



dVL 



dk 



I. 



it)akx{t)) (23) 



where k = k/|k| represents the direction of measurement. 
Since, according to Eq. (f22|) . the number of photons is 
proportional to pulse length as expected, it is more con- 
venient to calculate the number of photons scattered into 
dQ per unit of time 



d^N, 



ph 



dVldt 



(k,kL) 



n^Luld^ 

327r2A2eo?ic3 
d^K 



W(k)F(k-kL) 



dndt 



(k,kz 



one 
atom 



F(q) (24) 



where >V(k) = (l — (ed • Ck)^) is the dipole pattern of 
the emitted light and Cd is a unit vector in the direction 
of the dipole moment d that is determined by a polariza- 
tion of the probing laser. Eq. ([M)) shows that the angular 
distribution of the scattered light, apart from the contri- 
bution from the dipole pattern, is determined only by 
F(q). In addition, for a single atom F(q) = 1 and thus 
all the information about scattering from the system of 
N atoms is contained in F{q). Therefore, in the subse- 
quent sections, we can focus solely on the properties of 
F{q), keeping in mind that the remaining contribution 
is the same as for the scattering from a single atom. 



The generalization of our derivation to the case of in- 
teracting atoms can be performed in an analogy to the 
problem of neutron scattering from liquid helium [2^.[26j. 
If one applies the Born approximation, and eliminates 
adiabatically the excited state, one ends up with the re- 
sult identical to the one presented here. A similar ap- 
proach has been applied to the study of Raman scattering 
in the superradiant regime 27 1. 

Finally we note that our perturbative treatment ne- 
glects the effects of the momentum transfer resulting 
from the photon recoil in the process of light scattering. 
We assume, however, that the scattered light is weak 
and far-detuned, therefore we expect that the fraction of 
atoms which experience the photon recoil is sufficiently 
small, such that the atom statistics is not significantly 
affected. Moreover, in the presence of a tight trapping 
potential, such as a deep optical lattice, one finds that 
the scattering is recoilless |28| . which requires the trap 
size smaller than a wave length of the scattered light. 



III. SCATTERING FROM ULTRACOLD GAS 
OF BOSONS IN AN OPTICAL LATTICE 



The setup we consider is schematically plotted in 
Fig. [1] It consists of an ultracold gas of bosons con- 
fined in an optical cubic lattice of M sites. We assume 
a homogeneous system with an equal average number of 
atoms n = N/M in each site of a lattice. The periodic 
potential of the lattice reads 0]: 



Vp{x, y, z) ^ Vq (sin^ kpX 



kpZj 



(25) 



where kp is a wave vector of laser beams that are used 
to form the lattice and Vq is the potential depth. The 




FIG. 1: (Color online) Setup. An ultracold gas of bosons 
confined in an optical lattice is illuminated with a probing 
laser beam (yellow arrow) characterized by the wavevector 
k_L. The photons scattered into a selected direction (green 
arrow) are collected by a detector. 
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exact configuration of the probing beam and detectors 
will depend on a dimensionality of the lattice and will be 
discussed later. In order to use the results of the previ- 
ous section, we need to specify a single-particle basis. In 
the case of atoms confined in an optical lattice it is con- 
venient to choose the basis of Wannier functions Wm{r) 
that represent wave functions localized at single lattice 
sites m and are linear combinations of Bloch states. In 
our approach we consider only excitations within the low- 
est Bloch band, so in a limit of deep optical lattices the 
Wannier functions describe only the ground state wave 
functions in local potential wells. 



A. Deep lattice regime 

For a deep optical lattice, the Wannier states arc well 
localized within the sites of the lattice, and in equation 
we can restrict to optical transitions between the 
states localized at the same lattice sites: n = n' and 
m = m'. In this case F{q) simplifies to the following 
expression 

F(q) = E(n|e''i"-|n)(m|e-'i--|m) (^i^ffn^Lffm) 

n,m 

= l/o.o(q)|' e"^*'""-'"'"' (n„n™) (26) 



where 



ffiiffm and 



/„,™(q) = (n|e*'i--|m) = J d\ w*Jr)e'^''w^{r). (27) 

In analogy to the scattering of light into an optical cavity 
[2^ . we can define the classical part F'^^°-'^{q) and the 
quantum part i^'?'"^"*(q) of the function F(q) 



F^'-(q)^n2|/o^o(q)r 



(28) 



F'?"""*(q) = F(q)-i^"'"-'(q) 



(29) 



The former yields the classical amplitude of the scattered 
light Kcica)!^, whereas the latter represents the remain- 
ing quantum contribution that together with _F'^'°''(q) 

sum up to the total number of photons (cLl^xCikx^- We 

note that i^^''"'(q) has a form characteristic for a Bragg 
scattering and it is not affected by any statistical prop- 
erties of the ultracold gas of bosons. On the contrary, 
F''"°"*(q) is sensitive to the atom number statistics and 
thus enables us to investigate statistical properties of dif- 
ferent quantum states. 



IV. SCATTERING FROM A BOSE GAS IN 
ONE-DIMENSIONAL OPTICAL LATTICE AT 
ZERO TEMPERATURE 

The geometry of the system we investigate in this sec- 
tion is depicted in Fig. (5] We consider one-dimensional 
homogeneous optical lattice generated by two overlap- 
ping and counterpropagating laser beams characterized 
by the wavelength Ap. Atoms confined to the periodic 
potential are illuminated with a single laser beam with 
the wavelength Al- For the single particle basis that we 
have chosen the states 



(30) 



arc products of a Wannier function Wmiz) localized at 
lattice site m along z-direction, and a Gaussian func- 
tion ip^(x,y) in tightly confined, perpendicular direc- 
tion. For simplicity we assume the cylindrical symmetry 
^j_{x,y) = Tpj^{p). At zero temparature a gas of bosons 
in a periodic potential appears in two distinct quantum 
phases [2^ [SO] ■ When the tunneling process dominates 
over the on-site atom repulsion the system is found in 
the supcrfluid (SF) phase that is characterized by the 
presence of a global coherence and a non-zero order pa- 
rameter. In contrast, for the on-site interactions stronger 
than the tunneling rate, the system exhibits the Mott- 
insulator phase. In the latter case the global coherence 
is lost, while the on-site particle number is fixed and the 
on-site fluctuations are suppressed. 

For a Bose gas at zero temperature and deep in the 
MI regime, the on-site fluctuations and correlations van- 
ish: {n,n.nm') — {n-rn) ('^rn') = 0. Hcuce, the quantum 
part F'?"°"*(q) is zero identically, and the scattering is 
described by the standard Bragg pattern with charac- 
teristic set of maxima and minima, corresponding to the 
directions of constructive and destructive interference. In 




FIG. 2: (Color online) Setup. A quasi one-dimensional opti- 
cal lattice is illuminated with a probing laser set at an angle 
a t [— TT, 7r[. A detector is set at an angle (5 e [— vr, 7r[. 
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contrast, SF phase at T = exhibits nonzero fluctuations 

and correlations: {rimnm') - (rim) {nm') = n5mm.' - jf- 
Hence, apart from the similar behavior of the classical 
part F'^'"*(q) as for MI phase, the SF phase also gives rise 
to nonzero quantum component F''"''"*(q) which, within 
the deep lattice approximation (Eq. ([26]) ). is given by 
^g«ant(-q) ^ N \fo,o{<i)f ■ This offers a unique possibil- 
ity of a non-destructive measurement that allows one to 
distinguish between SF and MI phases p^ . 

Fig. [3] and Fig. [4] compare the scattering patterns from 
the SF and MI phases for the systems of M — 55 sites 
with different configurations of the probing laser and dif- 
ferent ratios of Ap to A^. At some characteristic angles 
corresponding to the Bragg scattering minima due to the 
destructive interference, the scattering from the MI state 
vanishes. In contrast, the scattering pattern from the 
SF state is nonzero at all angles, also in the directions 
where the classical component vanishes. We observe that 
a change of a ratio Ap/Ai affects the scattering pattern, 
in particular a number and positions of the highest peaks 
resulting from the constructive interference. 

We note that for large M the scattering pattern quickly 
oscillates and thus, in the realistic measurement, one 
would detect photons scattered in some finite solid an- 
gle dfl which is characteristic for the detector and that 
contains several interference fringes. Hence, we find it 



/J /J 
(a) a = 7r/4 (b) a = 7r/2 

FIG. 3: (Color online) Structure function i^(q) for SF (blue 
top curve) and MI (gray bottom curve) phases, in the deep 
lattice approximation, for a probing laser set at different an- 
gles ct, and a detector set at angle (5. Here, Vb = 15Jf?r, 
M = 55, = 3M, \p/\L = 1. The black line represents 
average distribution (|32|l . 




p li 
(a) Ap/Az, = 2 (b) Ap/Az, = 1/3 



FIG. 4: (Color online) Structure function F(q) for SF (blue 
top curve) and MI (gray bottom curve) phases, in the deep 
lattice approximation, for different ratios of Ap/A_L, and a 
detector set at angle (3. Here, Vo = IQEr, M = 55, iV = 3M, 
a = 7r/8. The black line represents the average distribution 

(ED). 



more appropriate to calculate the angular distribution of 
photons that is averaged over few neighboring maxima. 
The averaging does not affect the scattering pattern of 
SF phase, which is rather smooth, but it is important for 
MI phase. In ID optical lattice, the angular distribution 
of photons scattered from MI state is determined by 

i^^"(q) = l/o,o(q)rn^^^^4#^ (31) 
sm (^qdj 

where d denotes the translation vector of a ID lattice. 
Averaging over some finite solid angle around q contain- 
ing several maxima yields 



F^'{q) = |/o,o(q)|' ^ . 2n (32) 
2 sm (^qdj 

The above result is derived provided that the measure- 
ment is done not too close to the main maxima deter- 
mined by the directions of the constructive interference. 
As can be observed in Fig. [3] and Fig.|4l the averaged dis- 
tribution of the light scattered from MI phase still can 
be well distinguished from the scattering from the SF 
phase, and result for the averaged distribution re- 
mains approximately valid even close to the points of the 
destructive interference. 

While performing the deep lattice approximation in 
Eq. (pS)) we have dropped all the Franck-Condon fac- 
tors corresponding to transitions between states localized 
in different lattice sites. Obviously, with the decreas- 
ing lattice depth, the Wannier states begin to overlap 
between neighboring sites and we expect the nonlocal 
Franck-Condon factors to give larger contribution. In or- 
der to investigate this issue in Fig. [S] we compare the deep 
lattice approximation (pS)) with the result that includes 
summation over all pairs of the lattice sites in Eq. p8|) . 
For the clarity of presentation we show the results for 
a relatively small system oi M = 11 sites and a lattice 
depth V = lEr and V = 5Er for SF and MI phases, 
respectively, expressed in the units of the recoil energy 
Er = fi^k\/ {2m). We observe that even for the shallow 
lattice potential the nonlocal corrections give negligible 
contribution for the scattering from the SF state. More- 
over, for the MI state the nonlocal corrections are even 
smaller because of the deeper lattice required to achieve 
this phase. In Appendix VK\ we show that corrections due 
to the nearest neighbors in weak lattices are isotropic and 
scale as the total number of atoms N . Therefore, nonlo- 
cal corrections give rise to a scattering at the angles of 
destructive interference of the classical part. However, 
for typical lattice depths, the corresponding contribution 
is small and can be totally neglected for both quantum 
phases. 
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(a) SF scattering at Vq = lEr (b) MI scattering at Vq = 5Er 

FIG. 5; (Color online) Distribution of light scattered from SF 
(left panel) and MI (right panel) phases at zero temperature, 
versus the angle (3. Here, M = 11, N = 3M, a = 7r/2, and 
^p/^L ~ 1- The red curves present F{P) calculated in local 
approximation (|26p . whereas the blue curve (in the SF case) 
and the gray curve (in the MI case) show results including 
also non-local Franck-Condon factors as in Eq. (|18|) . In the 
case of scattering from the SF state at Vo = lEr one notices 
slight differences between the two curves in a vicinity of /3 = 
and /3 = ±7r. In the case of scattering from the MI state at 
Vo = ^Er the two curves are indistinguishable. 



STATISTICAL PROPERTIES AT FINITE 
TEMPERATURES 



annihilation and creation operators in momentum space 
1 



Ok 



1 



(36) 
(37) 



respectively, in which index m runs over all sites in a 
lattice. A period of the cubic lattice \s d ~ Xp/2 and a 
size of the system is equal to L = M^^^d. The periodic 
boundary conditions imply quantization of a wave vector: 
k = ^{nxTny^riz), where Uj are integer numbers ranging 
from - [M/2\ to [M/2\ [3l|. By rewriting Eq. ([33]) in 
terms of Uk and aj^., we find: 

7^ = ^ ekfllflk + ^ 4+k"«l<'-k"ak'ak (38) 



k,k'.k" 



where 



3 

ek = 6 J — 2 J ^ cos (kid) . 
=1 



(39) 



A. Bose-Hubbard model of an ultracold gas in a 
periodic potential 

As discussed in the previous section, the angular dis- 
tribution of the scattered light is determined by the occu- 
pation number statistics in lattice sites. We investigate 
the occupation number statistics within Bose-Hubbard 
(BH) model [1^, , considering only excitations within 
the lowest Bloch band. The BH Hamiltonian reads: 



B. Statistical properties of the superfluid phase 

The standard description of a weakly interacting Bose 
gas is based on the Bogoliubov approximation [32| and 
can be also applied for a superfluid phase in periodic 
potentials [33| . We perform the Bogoliubov approxima- 
tion to Hamiltonian (|38p . replacing the annihilation and 
creation operators in the zero quasi-momentum modes 
(k = 0) by C-numbers oq ~ aj sa V^o- By introducing 



n 



(m,m' 



the quasi-particle annihilation and creation operators 6k 
-f/^ rim (rim — 1) (33) ^^^^ respectively, which fulfill the standard bosonic 



where the first sum on the right-hand side is restricted 
to nearest neighbors only. The parameter 



d^rw*^{r) 



2m 



«^m'(r) (34) 



commutation rules 



(5k k' 5 a-nd are related to 



ak and aj^. by the canonical transformation 



"k Vk 



flk 

t 

k 



(40) 



we diagonalize Hamiltonian (j38[) obtaining 



is the hopping matrix element between neighboring sites 
m and m', and parameter: 



U= — I dV|w„,(r)r 



(35) 



corresponds to the strength of the on site repulsion of two 
atoms in a lattice site m. As before, Wm(r) is the single- 
particle Wannier's wavefunction of an atom occupying 
site m in a lattice. 

The BH Hamiltonian can be equivalently expressed in 
the momentum space, which is convenient for the analysis 
of the SF phase at finite temperatures and application 
of the Bogoliubov method. To this end we introduce 



n = Eo 



ftt^k^k^k- 

k 



(41) 



Here, Eq represents constant, ground-state energy term, 
hujk are the energies of the quasi-particle excitation spec- 
trum 



et + 2C/-f ek. 



m' 



(42) 



and real-valued coefficients Wk and Uk of the Bogoliubov 
transformation are given by 



1 /ek + C^ 

2 I fittik 



(43) 
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We note, that the excitation spectrum hujk depends on 
the condensate population Nq, which is known in the ht- 
erature as the BogoHubov-Popov spectrum, and is well 
suited to describe the statistics of a BEC at finite tem- 
peratures [1. [s^i- 

Below the critical temperature, when the condensate 
is macroscopically occupied, the occupation statistics of 
the quasi-particle modes is given by the Bose-Einstein 
distribution: 

(^k^k) =^^^^-^ ^ A, k^O, (44) 

{blb^blbi,, ) =/k/k' + 5k^k' {fi + /k) , k, kV 0, 

(45) 

where the value of the chemical potential fi is set to zero. 
This follows from the fact the condensate acts as a reser- 
voir of particles, and distributions of particles in excited 
modes are not restricted by the particle number conser- 
vation, which is consistent with the so-called Maxwell- 
Demon (MD) ensemble approximation [s^, [3^, [s^l • Ap- 
plying Bogoliubov transformation (PD|) . and Eqs. and 
(|45p we can easily find the mean occupation, fluctuations 
and correlations of the number of atoms in the quantized 
quasi-momentum modes: 

(nk) =(4 + «k)/k + t'L (46) 
{5W) = {ul + vlf {fl + /k) + ulvl (47) 
(nkTik') ^ulvl (1 + 4/k + 4/2) 4,_k' + («k) M (48) 

in which k 7^ k' 7^ and nk = ajj.ak is a particle number 
operator for a quasi-momentum mode k. Calculations 
of statistical quantities P5|) - P5)) within the Bogoliubov- 
Popov method require self-consistent determination of 
the mean condensate population A^o- First, iVo enters 
the excitation spectrum as a parameter. Second, it is 
determined by the statistics itself, 

No^N-Y, ("k) , (49) 

which yields 




(50) 

Finally, by transforming to the position space with the 
help of Eqs. ([55)1 and ([57]) . we evaluate single-site oc- 
cupation number statistics, i.e. single-site fluctuations, 
(^^rim) = (j^m) ~ "-^i ^'^d correlations, (nm^^m') — n'^, 
between each pair of sites in the lattice. At the final 
stage they are substituted into Eq. (pS]) determining the 
angular distribution of the scattered light. 

In Fig. [S] we present results for the statistics of the 
SF state realized in an optical lattice. We have cho- 
sen four different values of the strength of interactions 
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(a) Number of condensate atoms (b) Atoms number fluctuations 
in the lattice. in a single site of the lattice. 
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(e) Momemtum modes (f) Momemtum modes 

correlations. Here, k = (2,2,2). correlations. Here, k = (1,0,0) 

and k' = (1,1,1). 



FIG. 6: (Color online) Statistical properties of the SF phase 
in an optical lattice. Here, M — 11x11x11 and A'' = 3M. 
Different colors of the curves refer to different values of the 
parameter U: blue (solid) for U = 0, green (dashed) for U = 
1, yellow (dot-dashed) for [7 = 4, red (dotted) for [7 = 6. The 
selected values of U imply the following values of quantum 
depletion: 0%, 1.6%, 6.8%, 9.7%, respectively. The strength 
of interaction U and temperature T are expressed in units 
of J. 



U that correspond to quantum depiction ranging from 
to approximately 0.1, which should be proper in the 
regime of a weakly interacting gas where the Bogoliubov 
method is applicable. An upper limit of temperatures 
we consider is established by the conditions of validity 
of the MD ensemble approximation that, for sufRciently 
large systems, works well up to a temperature close to 
the critical temperature Tc- As expected, we observe 
that the on-site fluctuations increase monotonically up 
to Tc- In contrast, the correlations between popula- 
tions of different sites exhibit non-monotonic behavior 
that is strongly dependent on a distance between the 
considered sites. The phenomena can be understood by 
studying the behavior of these statistical quantities at 
small and at large temperatures. Readily, in the limit 
T ^ the on-site fluctuations and correlations follow 
the behavior presented in Section HVl ((J^rim) — ^ 
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and {nmn-m') — = —jq- On the contrary, at large 
temperatures the fluctuations and correlations can be 
described consistently within a model of N indistin- 
guishable particles distributed over M degenerate levels: 
(S^Um) = 'n? + n and (nm^im') — = 0. For both of the 
limits, the analytical expressions are derived under the 
assumption U / J — *■ 0, however the approximations work 
reasonably well also for finite values of the ratio U/J. 
The last two panels of Fig. [S] present correlations be- 
tween different modes in the momentum space. We note 
that the correlation between an excited and the conden- 
sate mode {n\^no) , exhibits a maximum at some moderate 
temperature. This follows simply from the competition 
between the process of thermal depletion of the conden- 
sate and a growth of the thermal fraction. Similarly, 
in case of correlations between two excited modes, we 
observe that at some temperature the initial growth of 
(nicTiic') is suppressed by decrease in population of these 
modes in favor of population of modes of some higher 
quasi- momentum. 

C. Statistical properties of the Mott-insulator 
phase 

We introduce grand canonical Bose-Hubbard Hamilto- 
nian /C: 

(m,m') m m 

(51) 

In order to describe quantum statistics of the Mott insu- 
lator phase at finite temperatures we adopt a mean-field 
decoupling approximation [H, . In analogy to the Bo- 
goliubov approach, we introduce a complex mean-field 
parameter ip = {gm) that can be physically interpreted 
as an order parameter that is nonzero if the system is 
superfiuid. Below the phase transition point, the sym- 
metry related to the gauge invariance of the phase is 
spontaneously broken and without losing generality we 
can assume that ip is real. The new parameter allows 
one to decouple the hopping term occurring in Eq. (|5ip 

ffLffm' = V' {gL + .9m') - V^^- (52) 

By performing this substitution we can decompose 
Hamiltonian ([5T|) into a sum of mean-field local Hamil- 

tonians /C^^: JC = ^/C^^, where 

m 

/C^^ = - 2DJ4> {g^ + gl) + 2DJiP'' - fxh^ 

+ \un^ {n^ - 1) . (53) 

At zero temperature, calculation of the ground-state en- 
ergy and its minimization as a function of the superfiuid 
order parameter ijj yields the phase diagram analytically 
[3^. However, for non-zero temperatures the model has 
no analytical solution, and one has to resort to numerical 



calculations. Namely, by diagonalization of Eq. (|53|) we 
calculate grand canonical partition function Z^-tp), 

Zii:)^Tr{e-^'^'^'}, (54) 

and on its grounds we determine the grand thermody- 
namic potential f2 (?/>), 

ni^) = -^\nZ{^). (55) 

Subsequently, by minimizing ^{ip) with respect to ip^ we 
obtain the equilibrium value of the order parameter that 
we use to calculate all the relevant thermodynamic quan- 
tities. In particular: 

(nM) =^Tr {n^e-^'c""} , (56) 
(S'nM) =^Ti-{nle-P'''^'}-n\ (57) 

In the homogeneous lattice (?T.m(/^)) ~ n, and this iden- 
tity is used to determine the value of the chemical poten- 
tial, for a given single site occupation n. 

Our mean-field approach assumes decoupling of dif- 
ferent sites, (rim^-m') — — (m 7^ m'), that agrees 
with the zero temperature statistics of the MI assumed 
in section IV. This, in general, is not valid at higher 
temperatures when the hopping between adjacent sites 
in non-negligible. Nevertheless, it is satisfied at smaller 
temperatures considered here. 

In Fig. [7] wc present the mean single-site occupation 
and fluctuations as a function of the chemical potential, 
and for different system temperatures. The calculations 
have been performed within the mean-fleld model for 
U/J = 128. One can see the existence of a character- 
istic temperature above which, the fiat steps in (nm(M)) 
disappear completely, and the curve becomes monoton- 
ically increasing. This crossover is accompanied by an 
appearance of nonzero fiuctuations for all values of fi 
presented in the plot. This corresponds to the transi- 
tion from the MI to the normal phase. In MI phase the 
system is infinitely compressible: d (n) /dfi = 0, while 
in the normal or SF phase the compressibility becomes 
finite: d (n) /dfi ^ 0. Since the on-site fiuctuations can 
be expressed as (S^n) = d (n) /d{Pij). therefore they can 
be nonzero only in the normal or SF phase. In order to 
distinguish between the normal and SF phases one can 
resort to the value of the order parameter ip. Finally, 
we note that our mean-field treatment neglects the ef- 
fects of correlations between different sites and the quan- 
tum fiuctuations. In the more accurate models that take 
these effects into account, the on-site fiuctuations become 
nonzero already in the MI regime, close to the boundaries 
with the SF or normal phases. 
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(a) Average atoms number in a (b) Atoms number fluctuations 
single site in the lattice. in a single site in the lattice. 

FIG. 7: (Color online) Statistical properties of Mott insulator 
phase for U = 128. Colors of the curves refer to different 
values of temperature being considered; blue (dotted) for T = 
0, dark green (dot-dashed) for T = 4, green (dashed) for 
T = G, orange (long-dashed) for T = 8 and red (solid) for 
T = 11. All the values of parameters are expressed in units 
of J. 



VI. FINITE TEMPERATURE SCATTERING IN 
THREE DIMENSIONS 

We consider three-dimensional cubic lattice and as- 
sume a sufficiently large value of the trapping potential 
depth Vo to neglect corrections from the nonlocal Franck- 
Condon factors. The setup is depicted in Fig. [51 The 
lasers creating an optical lattice (Ap, red arrows) are set 
along X, y and z axes. In general, the position of a prob- 
ing laser (A^, yellow arrow, characterized by angles 0^, 
Ol) and detector (green arrow, characterized by angles 
&d) can be optimized in order to minimize the con- 
tribution from the classical component in the vicinity of 
the direction of the measurement, cf. Appendix [Bl This 
is of particular importance in case of large lattices for 




FIG. 8: (Color online) Setup. A three-dimensional optical 
cubic lattice generated by lasers Xp (red arrows) is illuminated 
by a probing laser Al (yellow arrow) set at angles {4>l,0l)- 
A detector of scattered photons (green arrow) is ahgned in a 
direction {(j)d,Od)- 



which, due to a big number of interference fringes, the 
detector would collect the photons from several interfer- 
ence peaks. Here, though, we do not choose the optimal 
configuration, we consider some example geometry which 
not only sufficiently reduces an influence of the classical 
component but also offers relatively simple experimental 
realization. Namely, we choose the direction of the prob- 
ing beam along one of the diagonals of the lattice cube, 
ki = |ki|(^,— -^,^); and a detector centered around 

In analogy to the one-dimensional case, we expect 
that sharp differences in the intensity of light scattered 
from the SF and MI phases can be observed at angles 
for which the classical component F'^'"*(^,0) is negligi- 
ble. In Fig. [5] we present the zero-temperature structure 
function F(0, 0) for MI phase that is equivalent to the 
F^^"-'^ {(f) , 6) . Keeping in mind that the quantum compo- 
nent F^"°"*(0, 6) for the SF state is slowly varying and of 
order of N , we observe that {(pdjOd) is indeed a promis- 
ing direction for a measurement that can distinguish the 
two quantum phases. In Fig. [TUlwe corroborate this ob- 
servation by presenting scattering patterns for the super- 
fluid F^'^{(f),9) and Mott-insulator F'^"{(f>,9) phases at 
T = 0. The plots show cross-sections of F{(f), 6) along 
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FIG. 9: (Color online) Logarithm of the structure function 
F(q) for MI phase at zero temperature, as a function of spher- 
ical angles of detection, {9,(f)). The probing laser is set at 
{4>L,dL)- Bright regions correspond to directions in which 
a large number of photons is scattered. The yellow circle 
(pointed by the yellow arrow) indicates the direction of the 
probing laser (global maximum of number of scattered pho- 
tons). The green circle (pointed by the green arrow) refers 
to the direction of a detection k = 1^^^^ l("72 ' ""Tf ' which is 
discussed in details in the text. Here, M — 55, A'^ — 3M and 

Vo = 20Er. 
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(a) F(<t>,e = ea) 
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(b) F{<l> = <l>d,B) 



FIG. 10: (Color online) Zero-temperature F^^{e, 0) (top blue 
curves) and F^'\6,<f)) (bottom gray curves) with domains 
restricted to (f) (left figure) and 6 (right figure). The position of 
the photons detector [ifd, Od) is indicated with vertical dashed 
lines. Here, M = 55 x 55 x 55, iV = 3M, X^/Xl = 1 and 
Vb = 20Er. 



the planes of constant 6 and 0, respectively. Evidently, 
for the specific values of parameters we have chosen and 
for the assumed directions of the probing laser and of the 
detector, the difference of number of photons scattered 
from the SF and MI phases is of order TV ^ 10^ and thus 
should be readily measurable in experiment. 

We turn now to the thermal effects and their influence 
on the angular distribution of scattered photons. Ana- 
lyzing Fig. [11] and Fig. [1^ we observe the isotropic and 
monotonic growth of the intensity of scattered light with 
temperature, for both SF and MI phases. In the case of 
SF phase, this behavior can be explained on the grounds 
of the Eq. rewritten in the momentum representa- 
tion by means of transformation (|36p - (|37p . In particular, 
if we disregard anomalous averages while calculating ex- 
pectation values of the form (^o\^^(i\i2o\^,j^o,\s_4^ ^ i-c if we 
perform the approximation 

^a^^akaa^jflkj^ ~ '5ki,k2'5k3,k4 (ajkiakiaiL^aka 

+ '^kl,kl<5k2,k3 («ki«k2ak2«ki^ , (58) 

Eq. pH)) can be rewritten as: 



nq)-^l/o,o(q)|' 



^ (nk^k' 

k^^k' 



N{N - 1) 



E 

m 

g»(-k+k'+q)r, 



zqrn 



^l/o,o(q)r 



(59) 



When temperature increases, a number of particles oc- 
cupying excited modes grows (sec Fig.lH]), causing an in- 
crease in correlations terms {n\ji\^i). In consequence, the 
total intensity of the scattered light F^^[T) increases 
monotonically with temperature in any direction of mea- 
surement. We note that some correlations terms {n\!^n\^i) 
start to decrease above some characteristic temperature 
(see Fig. [S]). However this does not influence the to- 
tal structure function F^'^{T) that grows monotonically 
with T. 
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FIG. 11: (Color online) Structure function F((/), 6*) for SF 
phase of bosons in a three-dimensional optical lattice. The 
preferred position of the detector {(l>d,dd) is marked with 
vertical dashed lines. The plots show the cross sections 
along the constant 6 (left panels) and constant (right pan- 
els). The calculations have been performed for ®^Rb atoms, 
Xp = \l = 850nm, M = 11 x 11 x 11 and TV = 3M. The 
upper panels show results for constant Vb = 6.80i5r {U — 4) 
and increasing value of temperature: T = 0, 3, 5, 7, 9, 11 (or- 
dered from the bottommost to the topmost curve) with U 
and T being expressed in units of J. The bottom panels 
show results for constant temperature T — 10 and increasing 
value of trapping potential depth: Vb = S.GGiJr {U = 1), Vb = 
5.15Er {U = 2),Vo = 6.80Er ((7 = 4),^) = ((7 = 6), 

ordered from the bottommost to the topmost curve. The in- 
sets show the number of photons scattered in the direction of 
a detector {4>d,dd) versus temperature (top panel) and inter- 
action strength (bottom panel). 



Similarly, an increase in the interaction strength (an 
increase in the lattice potential depth, equivalently) re- 
sults in larger population of excited modes, partially due 
to increase in the quantum depletion. This behavior leads 
to a growth of the correlation terms in Eq. ([55]) and again 
to the monotonic increase in the full function F^^{U). 

In the case of MI phase, presented in Fig. [T^l the in- 
crease in the number of scattered photons F^^^(T) is fully 
determined by a temperature-driven growth of single-site 
fluctuations that have been presented in Fig. [T] 



VII. SUMMARY AND CONCLUSIONS 

We have investigated the scattering of a weak and far- 
detuned laser light from a system of ultracold bosons in 
an optical lattice. We have shown that the light scatter- 
ing can be used as a probe of the on-site quantum statis- 
tics, in particular fluctuations and correlations. Calcu- 
lating the statistics for the superfluid and Mott-insulator 
phases at finite temperatures, we have determined the 
angular distributions of the mean number of the scat- 
tered photons. The profiles of the scattered light are fully 
determined by the on-site particle number fluctuations 
and correlations and thus allow for an experimental ver- 
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FIG. 12: (Color online) Structure function F(0,6') for MI 
phase of bosons in a three-dimensional optical lattice. The 
preferred position of the detector {4>d,0d) is marked with ver- 
tical dashed lines. The plots show the cross sections along the 
constant 9 (left panels) and constant (j) (right panels). The 
calculations have been performed for *^Rb atoms, \p = \l = 
850nm, M = 55 x 55 x 55, iV = 3M, Vq = 18.3^^ {U = 128), 
/X = 320 and temperatures: T — 0,6,8, 10, 11 (ordered from 
the bottommost to the topmost curve), with U,T and ^ be- 
ing expressed in units of J. The inset shows the number of 
photons scattered in the direction of a detector {(pd, 6d) versus 
temperature. 



and SF phase, the coefficients /i.o(q) and /-i.o(q) rapidly 
tend to zero with the increasing lattice depth. 

Similarly, the nearest neighbor correction for super- 
fluid state reads 



AF^^(q) = (|/i,o(q)P + |/-i,o(q)P 
+ 2Re{/*o(q)/.i,o(q)} 



N - 



1 

1 

N 



(A2) 



A brief estimate leads to 

AF^^(q)«(|/i,o(q)P + |/-i,o(q)P) 



that, again, contains small coefficients /i,o(q) and 
/-i,o(q) rapidly decreasing with the lattice potential 
depth. 



ification of the present theoretical models describing the 
statistics in ultracold gases. For the 3D optical lattice we 
have determined the optimal geometry at which the con- 
tribution from the Bragg scattering pattern is minimized. 
We have shown that even at some non-optimal configura- 
tions, which can be more accessible from the experimen- 
tal point of view, this contribution is sufficiently small 
and allows one to measure the effects of quantum statis- 
tics. Our main conclusion is that by careful choice of the 
measurement geometry one can distinguish between dif- 
ferent phases, even at finite temperatures, and observe 
the effects of the finite temperature statistics of a Bose 
gas. 
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APPENDIX A: CORRECTIONS FOR WEAK 
LATTICE POTENTIALS DUE TO THE 
NONLOCAL FRANCK-CONDON 
COEFFICIENTS 

By applying the local approximation in the derivation 
of Eq. we have neglected contributions from the non- 
local Franck-Gondon coefficients. Here, we calculate the 
leading contribution from the neglected nearest-neighbor 
terms. 

In the case of the MI phase the expression is propor- 
tional to N: 

AF^"(q) = (|/i,o(q)P + |/-i,o(q)n (1 + n) N. (Al) 
Although it scales the same as the difference between MI 



APPENDIX B: OPTIMIZATION OF POSITIONS 
OF A PROBING LASER AND A DETECTOR IN 
THE 3D CASE 

In this appendix we derive the condition for optimal 
configuration of the probing light and of the photon de- 
tector, which lead to the minimal contribution from the 
classical amplitude of the scattered light. We start with 
the classical part of the structure function i^(q), defined 
in (121 



F^'-(q)=n2|/o,o(q)r 



(Bl) 



The label m, enumerates lattice sites: 
ymy + ZTTiz), with integer 



For a simple cubic 



lattice the summation can be easily performed 



F^'-(q)=n2|/o,o(q)|' J] 



sm 



i=x,y,z 



'\q,d) 



(B2) 



For the rest of the derivation we introduce a conve- 
nient parametrization of the vectors = kL{ax, a^, a^), 
k = kL{l3x,f3y,f3z), and q = kLirjx^riy, riz) describing the 
momenta of the incoming and scattered photons, and 
the momentum transfer, respectively. The dimensionless 
numbers a^, and /?i satisfy: \rji\ < 2, \ai\ < 1, and 
< 1 for i — x,y, z. Expressing the translation vector 
of the lattice and the wave vector of the laser in terms of 
the wavelengths: d = Xp/2 and — 2tt/Xl, we rewrite 
Eq. (|B2p in the following way 



F^'"^(q)=n'|/o,o(q)|' n 



sm 



(B3) 



Typically, the angular dependence of the Franck-Gondon 
factor |/o,o(q)|^ is rather weak, which follows from the 
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fact that the characteristic size of a single lattice site, 
given by a harmonic oscillator length associated with 
the potential well, is much smaller than the wavelength 
of the probing laser. In such conditions the scattering 
due to |/o,o(q)|^ is almost isotropic, and most of the 
angular dependence is determined by the interference 
term characteristic for the Bragg scattering. The func- 
tion sin^ (AIx) I sin^ (x) generating the interference pat- 
tern, takes the maxima at a; = titt, while the minimal 
amplitude of oscillations occurs in the middle between 
two neighboring maxima: .t = 7r(7i + i). In fact, the 
latter determines the desired condition for the measure- 
ment with the minimal contribution from the classical 
component: = 1 + 271^, where ni are integers, and 

i = x^y,z. For simplicity we further consider only the 
simplest case \p = Xl- Since \rji\ < 2, the only possibility 
is j]i = ±1, which leads to the following three conditions: 

(^3 - = ±1' for j = V' 2- (B4) 

The other two conditions are given by the conservation of 
the momenta of the scattered photons: |k| = |ki| — ki,, 
which results in 

la.P + + = 1/3x1' + + 1/3.1' = 1- (B5) 

By combining (jB4p and (jB5P , we obtain the following two 
equations determining the coordinates of k and k^ : 

|a.P + |%P + |a,|2 = l, (B6) 
|a,±lp + |ay±l|2 + |a, ±1|2 = 1. (B7) 



Readily, there are infinitely many solutions of the two 
above equations. All of them lie on a circle that is a 
common part of two spheres in the three-dimensional 
space. One of the possible solutions is given by the set 
of numbers: = fc^ (6 — v^, 6 — \/6, 6 + 2-\/6) /12 and 
k = fci (-6 - V6, -6 - Ve, 2^/6 - 6) /12. 

Finally, we note that for the optimal geometry deter- 
mined by Eqs. (jB6[) and (jB7[) . an average of the classical 
component F'^'°''(q) over a finite solid angle containing 
several interference peaks results in the three-dimensional 
analog of formula ([5^ 



F^"^^(q) = |/o,o(q)|'^. (B8) 



Here, we have applied the condition sin (^^id) = 1 that 
follows the conditions for the optimal choice of the mea- 
surement geometry. We stress that this result is derived 
for this particular geometry, and only in this case the 
sine squared factors average out to ^ independently in 
all three directions. 
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